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SUMMARY 

The guided electromagnetic waves that propagate along a lossless dielectric rod im- 
mersed in isotropic and uniaxial plasmas are determined. For the uniaxial plasma, only 
the case where the optic axis is alined parallel to the dielectric rod is considered. Both 
plasmas are described by the linearized momentum transport equation for a cold, colli- 
sionless, electron fluid of uniform density. Numerical results for the propagation con- 
stants of the waves are presented as a function of the rod and plasma parameters for both 
the circularly symmetric and dipole modes. Important aspects of the results are dis- 
cussed with particular emphasis on plasma diagnostic applications. 


INTRODUCTION 

During the past several years, interest in the measurement of plasma properties has 
been increasing. Specifically, investigations in energy conversion, electromagnetic 
propulsion, the reentry problem, and microwave amplifiers have created the need for 
accurately determining the parameters of a plasma. 

Several new methods for measuring electron density are currently being investi- 
gated. One promising new method consists of relating the plasma parameters to the 
measurable properties of electromagnetic waves guided by a dielectric rod in the plasma. 
Recently, Robson and Stewart (ref. 1) presented a theoretical expression for the electron 
density as a function of the phase shift of the dominant HE mode along a dielectric 
rod in a cold, collisionless, isotropic plasma. Their experimental measurements of the 
electron density with this technique were approximately 30 percent higher than those ob- 
tained with a double Langmuir probe. 



A more extensive investigation of wave propagation on a dielectric rod in a plasma 
was made by Medecki, et al. (refs. 2 and 3) for electron density measurements on the 
Apollo reentry flights. Their theoretical investigation considered a cold, isotropic 
plasma with collisions; however, numerical results were only presented for the TEq^ 
mode in the collisionless case. Experimental values were obtained by placing conducting 
end plates on a dielectric rod to make a cylindrical probe resonator from which the elec- 
tron density was obtained by observing the resonant frequency shift. Good agreement 
was reported between theoretical and experimental values. 

Kreuzer and Muller (ref. 4) also considered electromagnetic wave propagation on a 
dielectric rod in a plasma. Their theoretical and experimental investigation emphasized 
the HEj^ mode. They found that the phase shift of the electromagnetic waves can be 
made a weak function of the electron density if the rod radius is large and a strong func- 
tion if it is small. A probe was constructed from a rod of large radius by narrowing a 
short section to a small radius. The electromagnetic waves were guided into and out of 
the plasma with virtually no phase shift due to the plasma on the large portion of the 
probe, while in the narrowed region, a measurable phase shift due to the plasma oc- 
curred from which the localized electron density was determined. Their experimental 
values were in good agreement with Langmuir probe data. 

The previous investigators (refs. 1 to 4) have been concerned primarily with the 
TE Q1 or the HE ^ ^ mode. In this report, a formal analysis is presented for the complete 
set of modes that propagate along a dielectric rod in a plasma. Several of the lowest 
order circularly symmetric and dipole modes are investigated in detail to determine the 
modes and the ranges of the rod and plasma parameters that are most suitable for diag- 
nostic applications. 

Since a static magnetic field is present in many plasmas, the analysis is conducted 
for both isotropic and uniaxial plasmas to assess the effect of a magnetic field. In prac- 
tice, a uniaxial plasma occurs whenever the electron cyclotron frequency due to the static 
magnetic field is much greater than the source and plasma frequencies. Only the case 
where the static magnetic field is alined parallel to the dielectric rod is considered. 


MODEL 

The model (fig. 1) consists of an infinitely long, circular, cylindrical rod (of radius 
a and relative dielectric constant k) which is immersed in a cold, collisionless, homo- 
geneous plasma. In practice, the plasma will be somewhat inhomogeneous since the 
electron density will be depleted near the surface of the rod because of the formation of 
a plasma sheath. The sheath region, which can be approximated by a thin layer of free 
space, can easily be taken into account when k = 1. The net effect in this case is 
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Figure 1. - Dielectric rod in plasma. 


simply to increase the rod radius by the sheath thickness since the rod and sheath have 
the same value of dielectric constant. However, the effect of the sheath when k is other 
than 1 cannot be taken into account as readily. Thus, to simplify the analysis, the effect 
of the sheath will be omitted. 


GENERAL SOLUTION FOR FIELD 
Basic Equations 

The basic equations to be used are Maxwell's equations and the linearized momentum 
transport equation for a cold, collisionless, electron fluid with an applied static magnetic 
field. These equations are 


V x e = -jo>ju 0 h 

(1) 

V X h = qN Q v + jwe Q e 

(2) 

jwm e v = q(e + v x p Q H o ) 

(3) 


where 

e time -varying electric field 

H q static magnetic field 

h time -varying magnetic field 



ixi mass of electron 

e 

N Q equilibrium electron density 

q charge of electron 

v electron fluid velocity vector 

e Q permittivity of free space 

ju Q permeability of free space 

o> angular source frequency 

It has been assumed that the time variation of all nonstationary quantities is e^ wt , 
where co is sufficiently large that the ion motion can be neglected. Equations (1) to (3) 
can be combined to yield 

V x e = -ju>ju 0 h (4) 

Vxh = jcoe o K * e (5) 

If the static magnetic field is either weak so that a> c « co^ and u> c « co (isotropic 
plasma) or strong so that co » co and c o » co (uniaxial plasma), the relative per- 

mmm C P C 

mittivity tensor or dyad K for the static magnetic field in the z-direction is given by 

g=e M + ¥fl) + € pVz (6) 

where for an isotropic plasma and = 1 for a uniaxial plasma and 

6p relative dielectric constant of isotropic plasma, 1 - (co^/co^ 
co _ angular electron cyclotron frequency, Iqp H /m I 

C | U U vJ| 

2 2 / 
co square of angular plasma frequency, q N /e m 

unit vector in i m direction 

Note that equations (4) and_(5) also apply in an ordinary dielectric medium with a relative 
dielectric constant k if K is replaced by Id, where I is the unit dyad. 


General Solution of Basic Equations in Plasma 

The solutions sought in this problem are electromagnetic waves guided by the rod 
which vary as e~^ z along the z-axis, where /3 is the propagation constant. These so- 
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lutions for the electric and magnetic fields are derived by first finding the solutions for 
the z -component of the fields from which the transverse components are than obtained. 
The electric and magnetic fields can be expressed in component form as 


e - ( e t + e A' 


(7) 

Mv>w 


(8) 


where 

e time -varying axial electric field 

e^ time-varying transverse electric field 

h time-varying axial magnetic field 

z 

hj. time -varying transverse magnetic field 

Likewise, the operator V can be expressed as 




where V t is the transverse operator (in cylindrical coordinates V t = a r (9/3 r ) 

+ ( i d /r)(d/d 9 )). 

Substituting equations (7) and (8) into equations (4) and (5) and equating like compo- 
nents yield 


v t xH t = i w V P Vz 

(9) 

v t xh z £ z- j/3a z xS t = ja,e o e l®t 

(10) 

V t x i t = -ju>M 0 h z a z 

(11) 

VVz-^z'V-^t 

(12) 


Now, solving equation (12) for 1^ and substituting this result into equation (10) yield 



where 




( 14 ) 


k 

° c 


c = 


( £ o“o) 


1/2 


(speed of light in a vacuum) 


The transverse electric field e^. can be expressed in radial and azimuthal component 
form as 


4 /am 0h 3e > 

e r = -L[— 2 


a 


2 \ r 90 9r 


(15) 


and 


. / 9h o de \ 

e g = — (-WM o “ + ~ — ] 

9 a 2 \ ° 0r r 90 / 


(16) 


Next, substituting equation (13) into (12) to eliminate e^ yields 


e.\n rv 


a>M 0 « 


X e£ + /3a. x { am V, x h £ + /3£ x V . x e aj 

Z Z Z y O t ZZ Z t Z zy 


(17) 


The transverse magnetic field hj. can be written in radial and azimuthal component form 
as 


. /-we^e, 9e„ 9h . > 

h = J-f — * + /3 — * 

<* 2 \ r 90 9r ) 


(18) 


and 
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(19) 


h n = — |we„e 


9e 


o c l 




a 


dr 



The differential equation for h is obtained by substituting equation (13) into (11) to 
z 

eliminate e t ; this yields 



( 20 ) 


2 2 2 2 
where V? is the transverse Laplacian (in cylindrical geometry Vf = (9 /9r ) 

t 2 9 2 t 

+ (l/r)(3/9r) + (l/r )(9 /d0 )). Likewise, the differential equation for e is obtained by 
substituting equation (17) into (9) which yields 


ft 2 - - 0 

where 


6 = 



( 21 ) 


( 22 ) 


The general solutions for equations (20) and (21) are readily found to be 


\ = ^ A n K n (ar)ei n0 e jajt (23) 

n=-°o 

e z = X) < 24 > 

n=-oo 

where 

n azimuthal wave number 

IV 

K (x) modified Bessel function of n m order 

n ' 

A constant to be determined 

n 

B constant to be determined 

n 

Since the fields must vanish as r — the branch of the square root in equations (14) 
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and (22) that has a positive real part must be chosen for the radial wave numbers a and 
6. Substituting equations (23) and (24) into (15), (16), (18), and (19) yields the general 
solutions for the transverse components of e and h. 

General Solution of Basic Equations in Dielectric 

As previously indicated in the section Basic Equations, the fields in a dielectric can 
be expressed by the plasma equations if the relative dielectric tensor K is replaced by 
kf. Thus, equations (20) and (21) yield 

( v ? + “oX = 0 (25) 

and 

ft + “oX = ° (26) 

respectively, where the z -dependence of the fields is assumed to be of the form e~^ z 
and 



The general solutions to equations (25) and (26) are 

h z = X) C n J ,!<“o r > eln9ei " t < 28 > 

n=-<*> 

OO 

e - V DJ <orr)<J n< ’eJ“ t (29) 

z / -v n n o v 7 

n=-°o 


where 


J n (x) 

Bessel function of order 

C n 

constant to be determined 

D n 

constant to be determined 
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Either sign of the square root in equation (27) may be chosen. 

The transverse components for the fields in the dielectric can be obtained by substi- 

o 

tuting equations (28) and (29) into equations (15), (16), (18), and (19) with ot replaced 
2 

by ~ a 0 , and replaced by k. 


Dispersion Equation 

The boundary conditions require e z , h z , e^ , and h^ to be continuous at the surface 
of the rod. These conditions lead to the following dispersion equation 



J k< a o a > 

J n<“o a > , 


K^(aa) 


r 

W 


K k< 6a >" 


K n (aa) 

< 

k 


+ £i 

_ K n (5a) _ 

■> 


(30) 


where the prime refers to differentiation with respect to the total argument. 

The solutions for /3 are, in general, both real and complex; however, only the real 
values of /3 which satisfy the dispersion equation will be presented since this analysis is 
only concerned with unattenuated propagating waves. It should also be noted that equa- 
tion (30) is an even function of both /3 and n. Thus, the spectrum of waves that propa- 
gates along the positive z-axis is identical to that which propagates along the negative 
z-axis; likewise, the spectrum of left circularly polarized waves is identical to that of 
the right circularly polarized waves. Hence, it is not necessary to investigate explicitly 
the dispersion equation for negative values of /3 and n. 


NUMERICAL RESULTS AND DISCUSSION 

The results are presented starting with the simplest limiting case of a dielectric rod 
in free space. Next, the results for a dielectric rod in an isotropic plasma are pre- 
sented, and finally the results for a dielectric rod in a uniaxial plasma are considered. 

These results are presented in the form of curves for the dimensionless propagation 
constant /3a as a function of the dimensionless source frequency k Q a for various dimen- 
sionless plasma frequencies co p a/c. The values of the dielectric constant used are 1, 2, 
and 10. 
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Except in the circularly symmetric case (n = 0), the waves that correspond to the 
various solutions for /3a are hybrids; that is, all three components of the electric and 
magnetic field are present. These hybrid waves are denoted as either electric hybrid 
(HE) or magnetic hybrid (EH) . Each hybrid wave can be considered to be the sum of a 
transverse electric (TE) and transverse magnetic (TM) wave, with the TE component 
being dominant in the EH wave and the TM component being dominant in the HE wave. 

In order to determine whether a particular solution for /3a corresponds to an EH or 
HE wave, the limiting case where o^a/c is infinite is considered. This limit reduces 
the model to an ordinary, dielectric fiUed, circular, waveguide where all the waves are 
TE and TM. If a particular hybrid wave reduces to a TE in this limit, the hybrid wave 
is denoted by EH; if it reduces to a TM, it is denoted by HE. As previously mentioned, 
the circularly symmetric waves are not hybrids; all these solutions correspond to TE 
and TM waves. 

In addition, each wave is further denoted by the subscripts n and m which signify 
the azimuthal and radial orders of the wave, respectively. Values for these subscripts 
are selected so that in the limiting case where co^a/c is infinite they coincide with the 
subscripts of the corresponding circular waveguide modes. 


Dielectric Rod in Free Space 

The dispersion equation for a dielectric rod in free space can be obtained from equa- 
tion (30) with = €p = 1. For the circularly symmetric case (n = 0), the dispersion 
equation reduces to 


and 


where 


kQ? l aK o^l a ) wv* 

K 1 (a 1 a) J l( a o a ) 


(31) 


K^c^a) J i( Qf 0 a ) 


(32) 


“i a = [W 2 - M 2 ] 1/2 


(33) 
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It can be shown that solutions of equations (31) and (32) correspond to TM Qm and TE Qm 
modes, respectively. For the dipole case (n = 1), equation (30) reduces to 



The solutions of equation (34) correspond to both EH^ m and HE modes. 

Before the numerical results for the dispersion equations are presented, it is worth- 
while to consider the region in the (/3a - k Q a) plane where solutions can exist. For a di- 
electric rod in free space, solutions to the dispersion equations are possible only if the 
Kn( a i a ) function is real. This only occurs if the argument ot^a is real; consequently, 

/3a is bounded from below by the line /3a = k Q a. Moreover, it can be shown by consider- 
ing the sign of each term in the dispersion equations that /3a is bounded from above by the 
line /3a = k*^^k a. These bounds are shown in figure 2 and are valid for all values of n. 



Dimensionless source frequency, l^a 
Figure 2. - Region of solutions for dielectric rod in free space. 
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Figure 3. - Dielectric rod in free space. 


There are, of course, no solutions to the dispersion equations for a dielectric rod in 
free space if k = 1. Numerical results for k = 2 and k = 10 are presented in figure 3. 
Note that the TMq^ and the EH^ modes (as defined on p. 10) do not exist for the dielec- 
tric rod in free space. * 


■*ln the literature on surface waveguides (see, e. g. , ref. 5) the mode as de- 

fined in this report is often called the TM Q1 mode, the TM Q3 mode is often called the 
TMq 2 mode, etc. 
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The results show that all modes except the HE ^ ^ have a low frequency cutoff which 
occurs on the line /3a = k Q a. For the case n = 0, it can be shown that the only values for 
/3a and k Q a on this line which satisfy equations (31) and (32) are given by £* o a = |Wk o a^ 2 
oil / 2 ,, 

- ^ a > J = °om where a om is the m zero of V ff om^ (a 01 = 2 ‘ 405 ’ a 02 

= 5. 520 . . .)• Thus, the lower cutoff frequency for the TE Qm and ™ 0 ( m+ i) modes is 


k o a = 


om 


(k - 1) 


1/2 


(35) 


For the case n = 1, the lower cutoff frequency obtained from equation (34) for the 
EH l(m+l) and HE l(m + l) modes is 


k 0 a = 


g lm 

(k - 1)V2 


(36) 


where cr lm is the m th zero of ( a io = a n= 3.832. . .). There are no 

upper cutoff frequencies for modes on a dielectric rod in free space; as k Q a increases, 
the value of /3a for each mode asymptotically approaches k^^ 2 k 0 a. 


Dielectric Rod in Isotropic Plasma 


The dispersion equation for a dielectric rod in an isotropic plasma is given by equa- 
tion (30) with = £p. For n = 0 it yields 


ka 2 aK o (Q! 2 a ) , swv* 

E 1 ( Qf 2 a ) J 1 (%a) 


(37) 


for TML„, modes and 
om 


VW + VW 

Kj^a) Jj^a) 


( 38 ) 


for TE^ m modes where 
om 
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(39) 


a 2 a = [(/?a) 2 - e p (k 0 a) 2 ] 1/2 


For n = 1 it yields 


f-I 

a 2 BK 0 (a 2 a) * ^ 

+ 1 


\ 

r-v 

“ 2 aK 0 ( Q! 2 a ) ( x 

W) 2 

K^a) 

(v ) 2 

. J 1<V> 

f < 

l(v ) 2 

K 1 ( a 2 a ) 


+ k 

Q'aj(o'a) 

o o' 0 7 _ j 

ctf 

1 

1 . 1 

(v ) 2 

J i(« 0 a) 

-* 1 

WJ 

M M 2 


The region in the (/3a - k Q a) plane where solutions to the 


for EH lm and HE lm modes 
dispersion equations exist for an isotropic plasma is shown in figure 4. 

K 


It can be shown 


that if k Q a < co p a/c, all values of /3a are allowed 


k_a > o) p a/c, /3a is bounded from 

above by the line /3a = k*/ 2 k Q a and from below by the curve /3a = ^k Q a^ 2 - a/c^ ^ . 
These bounds are valid for all values of n. 



Figure 4. - Region of solutions for dielectric rod in plasma. 
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Numerical results for a dielectric rod in an isotropic plasma are presented in fig- 
ure 5. It should be noted from figures 5(d) to (i) that the presence of the plasma allows 
the TM Q1 and EH^ modes to propagate in contrast to the free-space case where they are 
cut off. The dispersion curves for these modes approach the line k Q a = 0 as o^a/c be- 
comes small. 

For low plasma frequencies, both the TMq ^ and EH^ modes are backward waves. 
This is particularly well illustrated in figures 5(d) to (f) for the TMq^ mode with 
o> a/c = 1. As the plasma frequency increases, these modes change from backward to 
forward waves. However, this transition does not occur at a discrete plasma frequency, 
but rather over a narrow range. If the dispersion curves for these modes are examined 
on a greatly magnified scale in this range, it is found that they are double valued with one 
value corresponding to a forward wave and the other to a backward wave. For the TMq^ 
mode, this narrow range can be approximated by a discrete plasma frequency which is 
given by 


c 


2.4 



(41) 


This equation was obtained from equation (37) by equating the value of k Q a for infinite 
/3a to the value of k Q a for /3a = k^'^a. For the case shown in figures 5(d) to (f) where 
k = 1, 2, and 10, these plasma frequencies are a^a/c = 2. 4, 2. 1, and 1. 8, respectively. 
No suitable criterion exists which yields a simple expression for approximating this nar- 
row range of plasma frequencies for the EHjj mode. 

As k Q a approaches the value 


ct> a 
P 


k 0 a = 


(k+ 1) 


1/2 


(42) 


the value of /3a for the TMqj and EH^ modes approaches infinity. This value of k Q a 
corresponds to a low-frequency cutoff for backward waves and a high-frequency cutoff for 
forward waves. Conversely, the high-frequency cutoff for backward waves and the low- 
frequency cutoff for forward waves cannot be expressed as simply since they involve the 
solutions of the transcendental dispersion equations (37) and (40). 

The dispersion curves for the remaining modes shown in figure 5 always correspond 
to forward waves. Except in the special case where k = 1, which is treated separately, 
these curves are similar for low values of plasma frequency such as co^a/c =1.0 since 
they are merely perturbations of the free-space case. For higher plasma frequencies 
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(a) Dielectric constant, 1. 



(d) Dielectric constant, 1. 


(b) Dielectric Constant, Z 


tc) Dielectric constant, 10. 



(e) Dielectric constant, Z (ft Dielectric constant 10. 


Figure 5. - Dielectric rod in isotropic plasma. 





(j) Dielectric constant, L 


Dimensionless source frequency, l^a 
(k) Dielectric constant, 2. 
Figure 5. - Concluded. 


(I) Dielectric constant, 10. 




such as o> p a/c = 10. 0, the curves approach those for a dielectric filled circular wave- 
guide; for co p a/c = °° the curves are, of course, identical to those for a dielectric filled 
circular waveguide. 

Each of these modes has a lower cutoff frequency but no upper cutoff frequency. As 
k Q a becomes large, the value of /3a for each mode asymptotically approaches k^^k a 
as in the free -space case. This is physically understandable since the effect of the plas- 
ma is negligible if k Q a » o; p a/c. For all these modes except the HE.^, the low- 
frequency cutoff first decreases with increasing aj p a/ c until a critical value for u> p a/ c 
is reached and then continuously increases with increasing o> p a/c. For values of o> p a/c 
below the critical plasma frequency, the low-frequency cutoff is given by 



for the TE Qm and TM Q ^ m+ ^ modes and 



(44) 


for the EH^ m+ ^ and HE^ m+ ^ modes. At the critical plasma frequency, cutoff occurs 
at k Q a = w p a/c, where co p a/c is given by 


to a 
P 

c 



(45) 


for the TE Qm and ™ 0 ( m+ i) modes and 

c " k l/2 


(46) 
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for the and EH i( m +l) modes> For values of o^a/c above the critical plasma 

frequency, the low-frequency cutoff for each mode is obtained by solving the appropriate 
transcendental dispersion equation with /3a = 0. For the HE mode, the low-frequency 
cutoff is zero for a^a/c = 0 and continuously increases with increasing plasma fre- 
quency. 

For the special case where k = 1, there are no solutions to the dispersion equation 
other than the TMq^ and EH n j modes unless the plasma frequency exceeds a critical 
value. It can be shown that this value is given by equations (45) and (46) for n = 0 and 1, 
respectively. For k = 1, all modes other than the TMq^ and EH n ^ modes have a lower 
but no upper cutoff frequency as in the case for general k. As k Q a becomes large, the 
value of /3a asymptotically approaches the value of k Q a. The lower cutoff frequency for 
each mode continuously increases as co^a/c increases, with the minimum value of the 
lower cutoff frequency being equal to the critical plasma frequency. 


Dielectric Rod in Uniaxial Plasma 

The dispersion equation for a dielectric rod in a uniaxial plasma with the optic axis 
(i. e. , direction of strong static magnetic field) alined parallel to the rod is given by 
equation (30) with = 1. For n = 0 it yields 

k(a' 1 a^ 2 K (6 1 a) or aJ (a a) 

V 1 l ° J_ + 0 0 0 = o (47) 

(G^K^a) J^a) 

for TM modes and 

om 

a.aK (a^a) a aj (a a) 

- 1 ° + 0 0 0 = 0 (48) 

K 1 (a 1 a) J i(“ 0 a ) 

for TE Qm modes where 

V = { € p [ <ftl)2 ■ (V)f /2 (49) 

and for n = 1 it yields 
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T“o aJ o<“o a ) 

“ 


ofjaKjo'ja) 




P»o a, o<“o*» 


“6 1 aK 0 (6 1 a) ' 

\ 

J 1<V> 

i 


K 1 (« 1 a) 


> < 

k 

J i<V> , 


K 1 (6 1 a) 

+ 1 









> 

(<v ) 2 



M 2 




(v ) 2 


(v ) 2 


i* 

- 


- 




- 


- 

-X 



for EHj m and HEj m modes. 

The regions in the (Pa - k Q a) plane where solutions to the dispersion equations can 
exist for the uniaxial plasma are shown in figure 6. It can be shown that if k Q a > V / c ’ 
solutions for all modes (i. e. , all values of n) exist in the region which is bounded from 
below by the line /3a = k Q a and from above by /3a = k^ // ^k Q a. These bounds are identical 
to those for the free-space case. When k Q a < cu^a/c only the circularly symmetric 
modes can exist. For the TE modes, the upper and lower bounds on /3a are 
Pa = k ' k Q a and /3a = k Q a, respectively. For the TM Qm modes, the upper and lower 



Figure 6. - Region of solutions for dielectric rod in uniaxial plasma. 
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bounds are 0a = k Q a and 0a = 0, respectively. It should be noted that if k = 1, only 

TM„ modes can exist and then only for k a < o> a/c. 
om o p 

Comparing equation (48) with (32) reveals that the dispersion equations tor the TE Qm 
modes are identical for the uniaxial and free -space cases. This result should be ex- 
pected since the electrons in a uniaxial plasma can only move parallel to the static mag- 
netic field which is in the z-direction, and the TE Qm modes, by definition, have no com- 
ponent of electric field in this direction. Thus, there is no interaction between the plas- 
ma and the TE^^ modes for the uniaxial case, 
om 

Numerical results for the TM Qm modes are given in figure 7. First, consider these 
results for k = 1. The TMqj mode for the uniaxial case exists for all values of co p a/c 
and has both a lower and an upper cutoff frequency as in the isotropic plasma case. The 
lower cutoff frequency can be obtained from equation (47) with 0a = 0; the upper cutoff 


frequency occurs along the line 0a = k Q a at k Q a = u) p a/c. The TM Q ( m+ ^ modes (m > 1) 
exist only if o) p a/c > Oj m . They also have both upper and lower cutoff frequencies. The 
lower cutoff frequencies for the TM Q ^ m+ ^ modes are given by equation (47) with 0a = 0, 
where the lowest possible value for the (m + 1)^ mode is k Q a = . The upper cutoff 

frequency occurs at k Q a = u> p a/c, where the value of 0a at the upper cutoff frequency is 

* - [M 2 - flm ) 2 ] 1/2 


Next, consider the case where the dielectric constant k is greater than 1; for this 
case, the TM om modes exist for all values of o) p a/c. Again, the TM Q1 mode has both 
an upper and lower cutoff frequency. The lower cutoff frequency is obtained as before by 


setting 0a = 0 in equation (47). The upper cutoff frequency occurs along the line 
0a = k Q a at k Q a = w p a/c if o> p a/c < a Q1 /(k - l) 1 ^ 2 , and at k Q a = °Qi/( k - l) 1 ^ 2 if 

o> p a/c > °Q-j/( k - l)*/ 2 . The results for the TM Q ^ m+ ^ modes (m > 1) differ from the 


TMqj modes since they have a lower but no upper cutoff frequency. Discussion of the 
lower cutoff frequencies can be divided into three cases. For o> p a/c < % m /(k - 1 ) 1/2 , 

the cutoff occurs on the line 0a = k Q a at k Q a = ° om /(k - l)^^ 2 . For 


CT om/ (k “ ^ 1/2 < ‘V' 70 < ff 
k Q a = a) a/c. Finally, for 
tion (47) by setting 0a = 0. 


l m /k* //2 , the cutoff occurs on the line 0a = k Q a at 
co p a/ c > CTj^/k^ 1 ^ 2 , the cutoff can be obtained from equa- 


The dispersion curves for the TM / m+ i\ modes (m > 1) are interesting since they 

1 /2 

have discontinuities or gaps in both 0a and k Q a. For °lm/ k < w p a /° 

< a o(m+l)/ (k ‘ 1)1//2 » a gap in 0a occurs at k Q a = u) p a/c. For a given u) p a/c in this 
range, the upper and lower values of 0a at the gap are the maximum and minimum 

values of |k^a) p a/c^ - MT 2 - co p a/c, respectively. Note that the gap vanishes 
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for OJ^a/c = cr lrn /(k - l) 1 / 2 . For o^a/c > °o( m+ i)/(k " 1) the gap occurs both in 
/3a and k Q a. The lower coordinates of the gap are /3a = a i( m+ i)/(k - l) 1 ^ and 

k Q a = a i( m+ i)/(k - 1) 1//2 ; the upper coordinates are /3a = 
and k Q a = o^a/c. 

Numerical results were also computed for the dipole modes in the uniaxial plasma. 
In the region of the (/3a - k Q a) plane where these solutions exist, the dipole solutions are 
practically indistinguishable from the corresponding dipole solutions for the free-space 
case. This is understandable since solutions for the case n = 1 cannot exist if 
k Q a < co^a/c (see fig. 6) where the strongest interaction between electromagnetic waves 
and the plasma normally occurs. If k Q a > co^a/c, the strong static magnetic field com- 
pletely eliminates the interaction of the plasma with the transverse electric field. Thus, 
only the axial component of the electric field can interact with the plasma. This inter- 
action is apparently small since the solutions are almost identical to those for the free- 
space case. 


[k(a> p a/c) 2 - (a lm ) 2 ] and 


APPLICATION OF RESULTS 


Dielectric rods in plasmas have been used to determine electron density in two basic 
ways. The first consists of terminating a short section of rod of length L with conduct- 
ing end plates to form a resonator (refs. 2 and 3). This structure will resonate when 
the value of k Q a is such that the corresponding value of /3a is an integral multiple of 
ira/L. These values of k Q a are indicative of the electron density. 

To facilitate the discussion of this method, it is useful to define a sensitivity S r as 
the rate of change of the resonant frequency with respect to the plasma frequency at a 
fixed propagation constant or . 


o = 9 (^) 
°r 




/3= constant 


9(k 0 a) 



/3a=constant 


(51) 


In order for the resonator method to be useful, the sensitivity must exceed a minimum 
value. For example, assume that a change in resonant frequency of 0. 1 percent can be 
resolved. Then, if a change in electron density of 5 percent is to be detected, the sensi- 

Q 

tivity S r must exceed 4x10" (co/Wp). 

The second method of determining the electron density consists of a direct measure- 
ment of the phase shift /3 l between points spaced a distance L along a dielectric rod in 
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a plasma (refs. 1 and 4). This measurement is normally conducted at a fixed source fre- 
quency k Q a. The sensitivity of this method is defined as the rate of change of the 
phase shift with respect to the plasma frequency at a fixed source frequency or 


S - 

<P 


3(<0p) 


_ L 3Q3a) 


w=constant 



k Q a= constant 


(52) 


To determine an approximate lower limit on the sensitivity S^, consider an example 
where the minimum phase shift that can be resolved is 0. 1 radian and where it is desir- 
able to detect a 5 percent change in electron density. This will require S to exceed 

4/ v 

The dynamic range for either of these methods is the range of o^a/c over which the 
sensitivity is sufficiently high to give the desired resolution in the electron density. 

The spatial resolution of the electron density measurement depends on the dimensions 
of the rod and the penetration of the electromagnetic field into the plasma. The spatial 
resolution in the axial direction is equal, of course, to the length L. The spatial resolu- 
tion in the radial direction depends on the radius of the rod a and the dimensionless ra- 
dial wave numbers aa and 6a. It can be shown from the asymptotic forms for the mod- 
ified Bessel functions in equations (29) and (30) that if &a > 1 and 6a > 1, the total ra- 
dial extent of the axial electric field is of the order of a[l + (l/6a)] and the total radial 
extent of the axial magnetic field is of the order of a[l + (1/aa)]. If cua « 1 and 
6a « 1, the radial extents of the axial electric and magnetic fields are generally quite 
large with respect to the radius of the rod. Thus, to obtain good spatial resolution it is 
necessary to select solutions where a a and 6a are at least 1 and a and L are small. 

The basic objective when applying the results, therefore, is to select the mode of 
propagation, values for the radius and dielectric constant of the rod, and values for 
either /3a or k Q a, that maximize the sensitivity and yet maintain the desired spatial 
resolution and dynamic range. 

First, consider the results for the isotropic plasma presented in figure 5. The 
TMqj and EHji modes appear to be very suitable for use in the resonator technique be- 
cause they yield a high sensitivity and a large dynamic range. Moreover, their sensi- 
tivities are nearly independent of the rod radius and plasma frequency if the value chosen 
for /3a is larger than the value of the largest dimensionless plasma frequency of inter- 
est. For large values of /3a, the sensitivity is the same for both the TMq^ and EHjj 

modes and can easily be computed using equation (51) giving S r = l/(k + 1)^^. This 
value is 25 times greater than the minimum value required for the example on page 23. 
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Since /3a is of the order of ira/L, both a and L can be small and still give the re- 
quired value for /3a. Thus, good spatial resolution can be obtained with these modes 
since the values of aa and 6a, which are of the order of /3a in this case, are large and 
a and L can be made small. 

The HE^ mode also appears to be useful in the resonator technique. For this mode, 
small values of /3a and o> p a/c give the largest sensitivity; however, the electromag- 
netic field may have a large radial extent since &a and 6a can be small if /3a and 
o)pa/c are small. For larger values of /3a, the sensitivity is generally very low unless 
o) p a/c is between 1 and 10. To obtain values in this range requires a prohibitively large 
rod radius for low plasma frequencies. For example, if the electron density N Q is 

10 15 per cubic meter, the rod radius must be of the order of 1. 7x10” * to 1. 7 meters; 
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however, if N Q is 10 per cubic meter, the rod radius must only be of the order of 
1. 7xl0 -5 to 1. 7xl0~ 2 meters. 

For the remaining modes, the value of /3a in most applications must be made larger 
than the value of /3a at the low-frequency cutoff for the free-space case to ensure that 
resonance will occur for all values of u> p a/c. Such a value for /3a will generally yield 
good spatial resolution but a low sensitivity unless cu p a/c can be made large. Again, 
this may result in a prohibitively large rod radius. Thus, the TM Q j and EH^ modes ap- 
pear to be the most suitable for use in the resonator technique - especially for the meas- 
urement of low plasma frequencies with good spacial resolution. 

All modes except the TMq^ and EH^ modes appear to be useful in the phase shift 
method. However, the source frequency for a particular mode must be selected suffi- 
ciently high so that the mode will propagate over the range of w p a/c of interest. Maxi- 
mum sensitivity generally occurs for values of o> p a/c in the range from 1 to 10. To 
achieve values in this range will require a large rod radius and hence give poor radial 
spatial resolution if the plasma frequency is low. If w a/c is restricted to be small, 

ir 

the sensitivity in the phase shift method can still be made arbitrarily large by increasing 
the length L (see eq. (52)). This can be done, however, only at the expense of de- 
creased axial spatial resolution. Thus, the phase shift method is generally not useful for 
measuring low electron densities with good spatial resolution. 

For the uniaxial case, only the TM Qm modes are useful; these modes are shown in 
figure 7. The TMq ^ mode appears to be usable in the resonance technique if the value 
chosen for /3a is less than the value of o> p a/c for the lowest plasma frequency of inter- 
est. In this region, the sensitivity is high, particularly for low plasma frequencies. 
Moreover, the sensitivity is nearly independent of the rod radius as in the isotropic case. 

The TMq 2 mode is usable in the resonator technique only for special cases since the 
dispersion curves are not a continuous function of co p a/c. The usefulness of the TMq 2 
mode is further restricted since the sensitivity is low unless the value of co p a/c is large 
which may require a large rod radius if the plasma frequency is low. Thus, the TM pi 


25 



mode appears to be the most suitable for use in the resonator technique for the uniaxial 
case. 

The uniaxial TMq ^ and TMq 2 modes appear to be very restricted in their usefulness 
in the phase shift method. The TMqj mode gives a high sensitivity and generally adequate 
spatial resolution but only over a very small dynamic range. It may be possible to extend 
the dynamic range, however, by operating at several different source frequencies. The 
TMq 2 mode gives a large dynamic range and adequate radial spatial resolution for low 
values of co p a/c; however, the sensitivity is very low unless the axial spatial resolution 
is sacrificed by making L large. For large values of o> p a/c, the sensitivity is gener- 
ally adequate, but the dynamic range is limited and the spatial resolution may be poor 
especially if the plasma frequency is low. Thus, the TMq ^ also appears to be the most 
suitable for use in the phase shift method for the uniaxial case. 

The parameter that has not yet been explicitly considered is k, the relative dielec- 
tric constant of the rod. A study of the results in figures 5 and 7 reveals that varying k 
produces two main effects. First, as k increases, the maximum sensitivity decreases 
for either technique discussed. This decrease is generally accompanied, however, by an 
increase in the dynamic range. Thus, the selection of the dielectric constant can be used 
advantageously to trade off sensitivity for dynamic range or the converse. The second 
main effect due to varying k is that either the source or resonant frequencies can be 
scaled with the scaling factor being approximately 1/k^^. This effect can be used in 
some cases to bring the frequency of interest within the range of available instrumenta- 
tion. Furthermore, by properly choosing k it may be possible to bring the radial extent 
of the field within acceptable limits. 

In addition to the selection of the rod and plasma parameters, several other factors 
must be considered before using electromagnetic wave propagation on a dielectric rod as 
a plasma diagnostic tool. One of these is mode identification. A general source for 
electromagnetic waves will excite not only the circularly symmetric and dipole modes, 
but also the quadrupole (n = 2) and higher modes. Thus, it is necessary to select a 
source configuration that will preferentially excite the mode or modes of interest so that 
they can be easily identified. This is particularly important when either the TMq^ or 
EH^ mode is used in an isotropic plasma since these modes along with the EH n ^ (n = 2, 

3, etc. , ) modes have their respective dispersion curves in the same region of the 
(/3a - k Q a) plane. 

Another factor to be considered is the effect of the inhomogeneous region in the plas- 
ma near the surface of the rod. As noted in the section MODEL, its effect for k = 1 is 
simply to increase the effective rod radius by approximately the thickness of this sheath. 
Since the sheath thickness is often unknown, it may be advantageous to operate where 
the relation between the wave and plasma parameters is not a strong function of the 
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rod radius. This is best illustrated by the resonator method with either the TMq j or 
EH^i mode. 


CONCLUDING REMARKS 

The guided electromagnetic waves that propagate along a lossless dielectric rod im- 
mersed in both an isotropic and uniaxial plasma were determined. Numerical solutions 
were presented for the propagation constants as a function of the parameters of the model 
for both the circularly symmetric and dipole modes. 

For the isotropic plasma, the TMq^ and EH^ modes were generally the most suit- 
able for use in the resonator method. They had good sensitivity, even for low electron 
densities, and a large dynamic range. For the phase shift method, all modes except the 
TM Q1 and EH ^ were acceptable; however, the values for the dimensionless plasma fre- 
quencies generally had to be in the range from 1 to 10 to give maximum sensitivity, which 
could result in poor spatial resolution if the plasma frequency is low. 

For the uniaxial case, only the TM Qm modes proved to be usable. The TMq ^ mode 
appeared to be the most suitable for use in the resonance technique if the value chosen for 
the propagation constant was less than the lowest value of interest for the dimensionless 
plasma frequency. The TM^ mode was only usable for special applications because of 
discontinuities in the dispersion curves. Both the TMq^ and TMgg modes were limited in 
their usefulness in the phase shift method. The TMq^ mode had a limited dynamic range 
which could, however, be expanded by operating at several different source frequencies. 
The TMq 2 mode generally had either inadequate sensitivity or poor spatial resolution un- 
less the electron density was high. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 16, 1968, 

125-24-03-05-22. 


27 



APPENDIX - SYMBOLS 


A, B 
n 7 n 

see eqs. (23) and (24) 

C, D 
n 7 n 

see eqs. (28) and (29) 

a 

radius of dielectric rod (see 
fig. 1) 

a i 

J-T_ 

unit vector in i ul direction 

c 

speed of light in vacuum 

e r 

time-varying radial electric 
field 

e z 

time-varying axial electric 
field 

e 0 

time-varying azimuthal elec- 
tric field 

e 

time -varying electric field 

®t 

time-varying transverse elec- 
tric field 

«o 

static magnetic field 

h r 

time-varying radial magnetic field 

h z 

time -varying axial magnetic 
field 

h 0 

time-varying azimuthal mag- 
netic field 

h 

time -varying magnetic field 

*t 

time -varying transverse mag- 
netic field 

I 

unit dyad 

J n ( x ) 

th 

Bessel function of n order 

K n (x) 

modified Bessel function of 
n^ 1 order 

K 

relative permittivity tensor of 


plasma (see eq. (6)) 

k relative dielectric constant rod 


k Q co/c (free-space wave number) 

L length of dielectric rod used as 

probe (for resonance tech- 
nique, L is distance between 
conducting end plates; for 
phase shift method, L is dis- 
tance over which phase shift is 
measured) 

in mass of electron 

e 

N q equilibrium electron density 

q charge of electron 

r radial coordinate (see fig. 1) 

S r sensitivity of resonator method 

(see eq. (51)) 

sensitivity of phase shift method 
(see eq. (52)) 

t time 

v electron fluid velocity vector in 

plasma 

z axial coordinate (see fig. 1) 

- (* 2 - ^lf 2 

“o « • ^f /2 

v [w 2 - (v ) 2 ] 172 

a 2 a jj/3a) 2 - £ p (V) 2 ] 1/2 

/3 propagation constant in z- 

direction 

• (v^ 2 - H*lf 2 

\ f P N 2 - M 2 ]} 1/2 
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I 


e Q permittivity of free space 

e„ relative dielectric constant of 

P 

isotropic plasma 

€j see eq. (6) 

P o permeability of free space 

9 coordinate in azimuthal direc- 

tion (see fig. 1) 

J.U 

(T m zero of J (x) 

om o' ' 


cr lm m th zero of Jjte) 

to angular source frequency 

(x> c angular electron cyclotron fre- 

quency 

cOp angular plasma frequency 

Subscripts: 

m order of mode in radial direction 

n azimuthal wave number 
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